Statistical regimes of random laser fluctuations 
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Statistical fluctuations of the light emitted from amplifying random media are studied theoreti- 
cally and numerically. The characteristic scales of the diffusive motion of light lead to Gaussian or 
power-law (Levy) distributed fluctuations depending on external control parameters. In the Levy 
regime, the output pulse is highly irregular leading to huge deviations from a mean-field descrip- 
tion. Monte Carlo simulations of a simplified model which includes the population of the medium, 
demonstrate the two statistical regimes and provide a comparison with dynamical rate equations. 
Different statistics of the fluctuations helps to explain recent experimental observations reported in 
the literature. 

PACS numbers: 05.40.-a,42.65.Sf,42.55.Px 



I. INTRODUCTION 



Optical transport in disordered dielectric materials can 
be described as a multiple scattering process in which 
light waves are randomly scattered by a large number of 
separate elements. To first approximation this gives rise 
to a diffusion process. The propagation of light waves 
inside disordered dielectric systems shows several analo- 
gies with electron transport in conducting solids [l| and 
the transport of cold atom gasses 0. A particularly in- 
teresting situation arises when optical gain is added to 
a random material. In such materials light is multiply 
scattered and also amplified. They can be realized, for 
instance, in the form of a suspension of micro particles 
with added laser dye or by grinding a laser crystal. Op- 
tical transport in such systems is described by a multiple 
scattering process with amplification. 

If the total gain becomes larger then the losses, fluc- 
tuations grow and these systems exhibit a lasing thresh- 
old. The simplest form of lasing in random systems is 
based on diffusive feedback [3| where a diffusion process 
traps the light long enough inside the sample to reach 
an overall gain larger then the losses. Interference effects 
do not play a role in this form of random lasing. Diffu- 
sive random lasing has been observed in various random 
systems with gain, including powders, laser dye suspen- 
sions, and organic systems [j, H, H, 0, Q ■ The behavior of 
such system shows several analogies with a regular laser, 
including gain narrowing, laser spiking and relaxation os- 
cillations 5, 9]. Reports in literature of complex emis- 
sion spectra from random lasers containing a collection of 
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spectrally narrow structures [HI UH 03 have triggered a 
debate about the possibility of lasing of Anderson local- 
ized modes in random systems [10( . Although Anderson 
localized modes can, in principle, form v ery interesting 
lasers resonators in a gain medium [ID, [lj], there has 
been no experimental evidence to date of random las- 
ing and localization in the same 3D random sample. In 
general, the observed spectra can be understood via a 
multiple scattering description based on the amplifica- 
tion of statistically rare long light paths that does not 
require localization or even interference [l5j . These emis- 
sion spectra exhibit a strongly chaotic behaviour, related 
to the statistical properties of the intensity above and 
around the laser threshold. 

Theoretical descriptions of light transport in amplify- 
ing disordered media and random lasing have been based 
so far on a diffusive mechanism 0, flrl_ 17|. using, for 
instance, a master-equation approach [18j |. To accom- 
modate the existence of localization related effects in 
the diffusion re gim e, 'anomalously localized states' have 
been proposed [IS |M III III iH IS • Other attempts 
to describe random lasing include Monte Carlo simula- 
tions [25|, [26j] , finite difference time domain calculations 
[27| . and an approach using random matrix theory [28| . 
A common feature of these studies is that the statisti- 
cal properties of a disordered optical system change with 
the addition of optical gain. Random lasers were found 
to exhibit, for instance, strong fluctuations of their laser 
threshold [IS [3(3] ■ It was also proposed that such sys- 
tems can exhibit Levy type statistics in the distribution 
of intensities [Hj] . 

In this paper we report on a detailed study of the sta- 
tistical fluctuations of the emitted light from random am- 
plifying media, using both general theoretical arguments 
and results from numerical studies. We find that the 
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characteristic scales of the diffusive motion of photons 
lead to Gaussian or power-law (Levy) distributed fluc- 
tuations depending on external parameters. The Levy 
regime is limited to a specific range of the gain length, 
and Gaussian statistics is recovered in the limit of both 
low and high gain. Monte Carlo simulations of a simpli- 
fied model which includes the medium's population, and 
parallel processing of a large number of random walk- 
ers, demonstrate the two statistical regimes and provide 
a comparison with dynamical rate equations. 

In Section |TT] we present some general arguments to 
explain the origin of the Levy statistics in random am- 
plifying media. In addition, we discuss the possibility of 
observing different statistical regimes. To check the va- 
lidity of the proposed general scenario, we define a simple 
stochastic model that is suitable for numerical simula- 
tions ( Section IIIip . The rate equation corresponding to 
its mean-field limit are introduced in Section HVl The 
results of Monte Carlo simulations are presented in Sec- 
tion [V] and discussed in the context of experimental re- 
sults in the concluding Section. 



II. STATISTICS OF THE EMITTED LIGHT 

Let us consider a sample of optically active material 
where photons can propagate and diffuse. Our descrip- 
tion is valid in the diffusive regime, hence we assume 
here that A is smaller then the mean free path I (A < £). 
The origin of the Levy statistics can be understood by 
means of the following reasoning. Spontaneously emitted 
photons are amplified within the active medium due to 
stimulated emission. Their emission energy is exponen- 
tially large in the path length I, i.e. 



1(1) = I exp{l/£ G ) 



(1) 



where we have introduced the gain length £q. On the 
other hand, the path length in a diffusing medium is a 
random variable with exponential probability distribu- 
tion 
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(2) 



where (I) is the average length of the photon path within 
the sample. The path length depends on both the ge- 
ometry of the sample and the diffusion constant D. A 
simple estimate of (I) can be provided by noting that 
for a diffusive process with diffusion coefficient D, (I) is 
proportional to the mean first-passage time yielding [32[ 



(3) 



where v is the speed of light in the medium and A is the 
smallest eigenvalue of the Laplacian in the active domain 
(with absorbing boundary conditions). For instance, A = 
q 2 with q — tt/L for an infinite slab or a sphere with L 
being the thickness or the radius, respectively 0. 



The combination of Eqs. (JTJ) and |J5J| immediately pro- 
vides that the probability distribution of the emitted in- 
tensity follows a power-law 
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Obviously the heavy-tail in Eq. ((4]) holds asymptotically 
and the distribution should be cut-off below some value 
Iq . The properties of the Levy distribution (more prop- 
erly termed Levy-stable) are well known [33j |. The most 
striking one is that for < a < 2 the average (I) ex- 
ists but the variance (and all higher-order moments) di- 
verges. This has important consequences on the statis- 
tics of experimental measurements, yielding highly irre- 
producible data with lack of self-averaging of samplc- 
to-sample fluctuations. On the contrary, for a > 2 the 
standard central-limit theorem holds, and fluctuations 
are Gaussian. 

The gain length Iq is basically controlled by popula- 
tion inversion of the active medium. Increasing the lat- 
ter, Iq and the exponent a (see Eq. (g])) decrease thus 
enhancing the fluctuations. At first glance, one may thus 
infer that the larger the pumping, the stronger the effect. 
On the other hand, £q is a time-dependent quantity that 
should be determined self-consistently from the dynam- 
ics. Indeed, above threshold the release of a huge number 
of photons may lead to such a sizeable depletion of the 
population inversion to force Iq to increase. It can then 
be argued that when the depletion is large enough, the 
Levy fluctuations may hardly be detectable. 

To put the above statements on a more quantitative 
ground we need to estimate the lifetime of the population 
as created by the pumping process. Following Ref. 0], 
we write the threshold condition as 



v/tc 



DA > 



(5) 



which is interpreted as "gain larger than losses" , the lat- 
ter being caused by the diffusive escape of light from the 
sample. Note that the condition r = along with Eqs. © 
and implies a = 1 at the laser threshold. 

For short pump pulses, the time necessary for the in- 
tensity to become large is of the order of the inverse of the 
growth rate r. When this time is smaller than the average 
residence time within the sample {l)/v, a sizeable ampli- 
fication occurs on average for each spontaneously emitted 
photon, leading to a strong depletion of the population 
inversion. In this case we expect a Gaussian regime where 
a mean field description is valid. The conditions for the 
Levy regime are therefore 1/r > {l)/v and a < 2 and can 
be written as: 



1 v . „ v 
2DA <iG< 2 DA 



(0) 



Note that the lower bound of the above inequalities cor- 
respond to a — 1/2. 

Without losing generality and for later convenience, let 
us focus on the case of a two-dimensional infinite slab of 
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thickness L. In Fig. [T] we graphically summarize equation 
© by drawing a diagram in the (L, Iq) plane. This rep- 
resentation allows to locate three different regions corre- 
sponding to different statistics. For convenience, the line 
corresponding to the threshold a = 1 is also drawn. The 
three regions of statistical interest displayed in Fig. [T] are: 

Subthreshold Levy: weak emission with Levy statistics 
with 1 < a < 2 (shaded region in Fig. [1] above the laser 
threshold line); 

Suprathreshold Levy: strong emission with Levy statis- 
tics with 1/2 < a < 1 (shaded region in Fig. [1] below the 
laser threshold line); 

Gaussian: a < 1/2 strong emission with Gaussian 
statistics, a > 2 weak emission with Gaussian fluctua- 
tions. 

Note that the first region corresponding to a nonlasing 
state, should display anomalous fluctuations as "precur- 
sors" of the transition. It should be also emphasized 
that the boundary between Levy and Gaussian statistics 
is not expected to correspond to a sharp transition (as 
displayed in Fig. [1]) but rather to a crossover region. 
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FIG. 1: (Color online) Different statistical regimes of fluc- 
tuations of a random laser with a two-dimensional slab ge- 
ometry of thickness L. For comparison with the simulations 
reported below, all quantities are expressed in dimensionless 
units I = 1, v = 1. The symbols correspond to the parameters 
of Fig. [2] 



III. NUMERICAL MODEL 

In order to provide evidence of the theoretical consider- 
ations presented in Section[nJ we introduce a general, yet 
easy to simulate, model of random lasing. We consider a 
sample partitioned in cells of linear size I. Specifically, we 
deal with a portion of a two-dimensional square lattice. 
Thus the center of each cell is identified by the vector in- 
dex r = (x, y), with x, y integers. In the following we will 
deal with a sample with a slab geometry i.e. 1 < x < L, 
1 < U < RL. The total number of lattice sites is thus RL 2 
where R defines the slab aspect ratio. Periodic boundary 
conditions in the y direction are assumed. 



Within each cell we have the population N(r, t) of exci- 
tations. We consider an hypothetical three-level system 
with fast decay from the lowest laser level. If the pop- 
ulation in the latter can be neglected we can identify N 
as the number of atoms in the excited state of the lasing 
transition. 

Isotropic diffusion of light is modeled as a standard 
random walk along the lattice sites. The natural time 
unit of the dynamics is thus given by At = Ijv. We 
choose to describe the diffusion dynamics in terms of a 
set of M walkers each carrying a given number of pho- 
tons ni, . . . tim- This may be visualized as an ensemble 
of "beams" propagating independently throughout the 
sample. Each of their intensities changes in time due to 
processes of stimulated and spontaneous emission. A ba- 
sic description of those phenomena can be given in terms 
of a suitable master equation [lU, HH that would require 
to take into account the discrete nature of the variables. 
To further simplify the model we consider that the popu- 
lation and number of photons within each cell are so large 
for the evolution to be well approximated by the deter- 
ministic equations for their averages. In other words, the 
rate of radiative processes is much larger than that of 
the diffusive ones and a huge number of emissions occurs 
within the time scale At [35| . With these simplifications 
N and n can be treated as continuous variables. Alto- 
gether the model is formulated by the following discrete- 
time dynamics: 

Step 0: pumping - The active medium is excited ho- 
mogeneously at the initial time i.e. N (r, 0) = Nq. The 
value Nq represents the pumping level due to some ex- 
ternal source. The initial number of walkers is set to 
M = 0. 

Step 1: spontaneous emission - At each time step and 
for every lattice site a spontaneous emission event ran- 
domly occurs with probability jNAt, where 7 denotes 
the spontaneous emission rate of the single atom. The 
local population is decreased by one: 



N 



N - 1 



(7) 



and a new walker is started from the corresponding site 
with initial photon number n = 1. The number of walkers 
M is increased by one accordingly. 

Step 2: diffusion - Parallel and asynchronous up- 
date of the photons' positions is performed. Each walker 
moves with equal probability to one of its four nearest 
neighbours. If the boundaries x — 1, L of the system are 
reached, the walker is emitted and its photon number 
n out recorded in the output. The walker is then removed 
from the simulation and M is diminished by one. 

Step 3: stimulated emission - At each step, the photon 
numbers ni of each walker and the population are up- 
dated deterministically according to the following rules: 



N 



(1+jAtN) n t 
(1 - ~/At nf) N 



(8) 



where N is the population at the lattice site on which 
the i— th walker resides. 



4 



Stochasticity is thus introduced in the model by both 
the randomness of spontaneous emission events (Step 1) 
and the diffusive process (Step 2). Note that the model 
in the above formulation does not include non-radiative 
decay mechanisms of the population. Furthermore, no 
dependence on the wavelength is, at present, accounted 
for; in general 7 = 7(A). 

The initialization described in Step is a crude mod- 
elling of the pulsed pumping employed experimentally. 
It amounts to considering an infinitely short excitation 
during which the samples absorbs N photons from the 
pump beam. As a further simplification we also assumed 
that the excitation is homogeneous throughout the whole 
sample. More realistic pumping mechanisms can be eas- 
ily included in this type of modeling. More importantly, 
as we are going to study the time dependence of the emis- 
sion, this type of scheme applies to the case in which the 
time separation between subsequent pump pulses is much 
larger than the duration of the emitted pulse (i.e. no re- 
pumping effects are present). 

Steps 1-3 are repeated up to a preassigned maximum 
number of iterations. The sum of all the photon num- 
bers of walkers flowing out of the medium at each time 
step is recorded. The resulting time series is binned on 
a time window of duration Tw to reconstruct the output 
pulse as it would be measured by an external photon- 
counter. This insures that each point is a sum over a 
large number of events and makes the comparison with 
ensemble-averaged results of the preceding Section sensi- 
ble. 

It should be emphasized that, although each walker 
evolves independently from all the others, they all in- 
teract with the same population distribution, which, in 
turn, determines the photon number distributions. In 
spite of its simplicity, the model therore describes these 
two quantities in a self-consistent way. 

For convenience, we chose to work henceforth in dimen- 
sionless units such that v = 1, I = 1 (and thus At = 1). 
The only independent parameters are then 7, the initial 
population Nq (i.e. the pumping level) and the slab sizes 
L, RL. 



1/4 in our case. 

For simplicity, let us consider the case of a laterally 
infinite slab (R —> 00) in which both TV and / depend on 
the y coordinate, only. Absorbing boundary conditions 
are imposed, 1(0, t) = I(L + l,t) = 0. The integration is 
started from the same initial conditions of the stochastic 
simulations, namely N(x, 0) = No, I(x, 0) = 0. 

As a first remark, we note that the threshold condition 
© applies to (fTT)|) upon identifying 



q = 



L + l 



(11) 



We can thus define a critical value of the initial popu- 
lation N c — Dq 2 /-f. For No < N c the total emission is 
very low being due to spontaneous processes that are only 
weakly amplified. On the contrary, for No > N c strong 
amplification occurs: the number of photons within the 
sample increases exponentially in time at a rate given by 
Eq. ©, r = jN - Dq 2 = 7(^0 - N c ). After the pulse 
has reached a maximum and the population is depleted, 
the emission decreases strongly. An estimate of the de- 
cay time of the pulse is given by solving the linearized 
equations around the stationary state N = 0, / = 0. 
A straightforward calculation yields that the long-time 
evolution is approximated by N(x,t) = N q (t)sxa.(qx), 
I(x,t) = I q (t) sin(gx) where 



N q (t) = Aexp(-Tt) 

7 



(12) 



I q {t) = A 



Dq 2 — 7 



exp(-7i) + B exp(-Dq J t), 



with A, B being suitable time-independent amplitudes. 

The above results have been checked by comparing 
them with the numerical solution of (|10p obtained by 
simple integration methods of ordinary differential equa- 
tions. In particular, we checked that both the rise and 
fall rates of the emission pulses (see the figures in the 
next section) are consistent with the expected values of 
r and Eqs. (fTB")) , respectively. 



IV. MEAN-FIELD EQUATIONS 

Before discussing the simulation of the stochastic 
model it is convenient to present some results on its 
mean- field limit. When both the population and photon 
number are large we expect the dynamics to be described 
by the rate equations for the macroscopic averages. This 
means that, up to relatively small fluctuations, the indi- 
vidual realization of the stochastic process should follow 
the solution of @, [3 



N = -jN(I+l) 
I = DAI + ~/N(I 



1) 



(9) 
(10) 



where 7(r, t) is the number of photons in each cell, A 
denotes the two-dimensional discrete Laplacian and D = 



V. MONTE CARLO SIMULATIONS 

In this Section we report the results of the simulation of 
the stochastic model. Preliminary runs were performed 
to check that lasing thresholds exist upon increasing of 
either the pumping parameter Nq and the slab width L. 
The values are in agreement with the theoretical analysis 
presented above. In addition, checks of relations ([2]) and 
© have been performed. 

As explained in Section ITTT1 we monitored the outcom- 
ing flux (per unit length) <j> as function of time. The re- 
sults are compared with the corresponding quantity eval- 
uated from Eqs. ([TO]) . In this case, <j> is defined from the 
discrete continuity equation to be 



D 
~2 



I(l,t) + I(L,t) 



(13) 
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FIG. 2: (Color online) The photon flux (per unit length) as 
a function of time for a sing le shot, for No = 2 x 10 9 (a); 
No = 5 X 10 9 (b) and No = 50 x 10 9 (c). Smooth red lines 
are the mean-field results, evaluated inserting the solutions 
of Eqs. (|10[) into Eq. I|13JI. For both curves, data have been 
binned over consecutive time windows of duration Tw = 10. 
Note the difference in the vertical-axis scales. 



The factor 2 comes from taking into account the contri- 
bution from the two boundaries x = 0, L of the lattice. 

The results of Monte Carlo simulation for a lattice with 
L = 30, R = 20 (18000 sites) and 7 = 10~ 12 (yielding 
N c = 2.5673 x 10 9 ) are reported in Fig.0 The three cho- 
sen values of Nq are representative of the three relevant 
statistical regions depicted in Fig. [TJ they correspond to 



£ G = 500 (Subthreshold Levy), l G = 200 (Suprathresh- 
old Levy) and £ G — 20 (Gaussian), respectively (see the 
triangles in Fig. [TJ. In the first two cases, the total emis- 
sion is highly irregular with huge deviations from the 
expected mean-field behavior. Above the lasing thresh- 
old (Fig. [Tb) single events ("lucky photons") may carry 
values of rij up to 10 10 . The resulting time-series are 
quite sensitive to initialization of the random number 
generator used in the simulation. On the contrary, in the 
Gaussian case (Fig. [5b) the pulse is pretty smooth and 
reproducible, except perhaps for its tails that, however, 
have a much smaller relative intensity. 

The evolution of the population N displays similar fea- 
tures. We have chosen to monitor the volume-averaged 
population 

^E^*) ( 14 ) 

r 

normalized to its initial value for a better comparison. 
Fig. [3] shows the corresponding time-series for the same 
runs of Fig. O Again, large deviations from mean-field 
appear for the first two values of No. The inset shows 
that in correspondence with large-amplitude events the 
population abruptly decreases yielding a distinctive step- 
wise decay. 

The non-smooth time decay is accompanied by irregu- 
lar evolution in space. Indeed, a snapshot of N(r,t) re- 
veals a highly inhomogeneous profile (see Fig. [4]) . Light 
regions are traces of high-energy events that locally de- 
plete the population before exiting the sample. 

For the Gaussian case, Fig. [5b, similar considerations 
as those made for the corresponding pulse apply. Note 
that now the population level decays extremely fast. It 
reaches 10% of its initial value at t » 600 which is only 
twice the average residence time within the sample. This 
means that photons emitted after a few hundreds time 
steps have hardly any chance to be significantly amplified 
(i.e. £ G has become too large). 

To check directly the validity of the power-law distri- 
bution ^ we computed the histograms of the photon 
number n ou t for each and every collected event during 
the whole simulation run (i.e. the same time range as 
in Fig. [2]). The result are given in Fig. [5] for three val- 
ues of N for which the Levy distribution ^ is expected 
to occur. A clear power-law tail extending over several 
decades is observed. Note that the middle curve cor- 
respond to the threshold value No = N c for which we 
expect a = 1. Remarkably, the values of the exponents 
measured by fitting the data are in excellent agreement 
with the definition of a (see inset of Fig. [5]) • As predicted, 
no meaningful value smaller than a = 1/2 is obtained 
from the data. It should be noted that we are deal- 
ing with a non-stationary process and the results may 
thus, in principle, depend on the observation time. To 
check for a possible time-dependence of the statistics, we 
considered a four-times longer simulation for the case of 
Fig. Wp an d divide the resulting time series in four consec- 
utive parts. Each of the resulting histograms are almost 
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time 

FIG. 3: (Color online) The normalized volume-averaged 
atomic population as a function of time for a single shot 
and for the same values of No as in Figs. [2^,-c (upper to 
lower curves respectively). Solid lines are the mean-field re- 
sults evaluated inserting the solutions of Eqs. (|10[) into (|14l) . 
The inset shows a magnification of the middle curves (case 
N = 5 X 10 9 of Fig. [2b). 




FIG. 4: A gray-scale plot of the atomic population distribu- 
tion along a portion of the lattice for t — 10 5 , No = 5 x 10 9 . 
White regions correspond to small values of TV. 



indistinguishable confirming that the underlying process 
is almost stationary at least on this time scale. 

Finally, to further elucidate the differences between the 
two types of statistical regimes, we performed a series of 
simulations increasing the number of lattice sites. For 
comparison, we kept L = 30 fixed and increased the as- 
pect ratio R up to a factor 4. In this way, we increased the 
number of walkers accordingly. For the Gaussian case, we 
did observe the expected reduction of fluctuations around 
the mean-field solution. On the contrary, the wild fluc- 
tuations of the Levy case were hardly affected. This is a 
further confirmation of the scenario discussed in Section 

VI. DISCUSSION 

Based on heuristic arguments, we have shown in Sec- 
tion |TT] that, depending on the value of the dimensionlcss 
parameter DA£g/v, the fluctuations in the emission of a 
random laser subject to short pump pulses can be dras- 
tically different. In a parameter region extending both 
above and below threshold, the intensity fluctuations fol- 
low a Levy distribution thus displaying wild fluctuations 
and huge differences in the emission from pulse to pulse. 




photon 



FIG. 5: (Color online) Histogram p(n ut) of the emitted pho- 
ton numbers n ou t for the two values of No of Figs. [2^, and b 
(lower and upper curves) and No = N c (middle) . This middle 
curve has been vertically shifted for clarity by roughly a factor 
5. The inset reports the values of a obtained by a power-law 
fit of the histograms as a function of the gain length £g as de- 
fined by Eq. {TTJ. The thin solid line is the theoretical curve 
as computed by Eqs. (|3]) and Q. The dashed vertical line 
marks the lasing threshold. 



In the suprathreshold case, such features have been in- 
deed observed in experiments [3l|. Some highly irrepro- 
ducible emission with lack of self-averaging and very ir- 
regular behaviour has been also detected [361 ] - 

The exponent a of the Levy distribution can be tuned 
upon changing the pumping level but it must be somehow 
bounded from below (a > 1/2) as a further crossover 
to a Gaussian statistics is attained. Indeed, far above 
threshold, when the gain length is very small, a large and 
fast depletion of population occurs (saturation). This 
hinders the possibility of huge amplification of individual 
events. In this case all photons behave in a statistically 
similar way. As a consequence, the statistics is Gaussian 
and a mean-field description applies again. 

The above considerations have been substantiated by 
comparison with a simple stochastic model. It includes 
population dynamics in a self-consistent manner. In the 
Levy regions, the simulation data strongly depart from 
the predictions of the mean-field approximation due to 
the overwhelming role of individual rare events. As a con- 
sequence, the evolution of the population displays abrupt 
changes in time and is highly inhomogeneous in space. 

To conclude this general discussion we remark that the 
width of the Levy region as defined by inequalities 
and depicted in Fig. Q] is of order L 2 . Since in our simple 
model, £q is inversely proportional to the pump parame- 
ter (see Eq. (TTTT) ). the interval of No values for which the 
Levy fluctuations occur shrinks as 1/L 2 . Therefore, the 
larger the lattice, the closer to threshold one must be to 
observe them. 

The existence of different statistical regimes, their 
crossovers and their dependence on various external pa- 
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rameters enriches the possible experimental scenarios. 
The emission statistics of random amplifying media has 
diverging moments in a finite region of parameters ex- 
tending across the threshold curve. Our theoretical work 
has shown that, depending on size, geometry, pumping 
protocols etc., the emission of random lasers may change 
considerably This general conclusion should be a useful 
guidance in understanding past and future experiments 
on random amplifying media. 
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